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Summary 

Coefficients of viscosity of various types of gas 
mixtures, including simulated natural-gas samples, 
have been measured at atmospheric pressure and 
room temperature using a modified capillary tube 
method. Pressure differences across the “straight” 
capillary tube section of a thermal mass flowmeter 
were measured for small, well-defined, volume flow 
rates for the test gases and for standard air. In this 
configuration, the flowmeter provides the volumet- 
ric flow rates as well as a well-characterized capillary 
section for differential pressure measurements across 
it. The coefficients of viscosity of the test gases were 
calculated using the reported value of 185.6 pP for 
the viscosity of dry air at 300 K. The coefficients 
of viscosity for the test mixtures were also calculated 
using Wilke’s approximation of the Chapman-Enskog 
(C-E) theory. The experimental and calculated val- 
ues for binary mixtures are in agreement within the 
reported accuracy of Wilke’s approximation of the 
C-E theory. However, the agreement for multicom- 
ponent mixtures is less satisfactory, presumably be- 
cause of the limitations of Wilke’s approximation of 
the classical dilute-gas state model. 

Introduction 

Measurement of viscosity of pure gases and gas 
mixtures at atmospheric pressure is of great inter- 
est in fluid flow and heat transfer studies. Such 
information is not always readily available for the 
gas mixtures, particularly for hydrocarbon mixtures. 
In recent studies (ref. 1) of heats of combustion 
of natural gas and mixtures of gaseous hydrocar- 
bons, reliable information about their viscosities was 
needed. In the absence of reliable and accurate ex- 
perimental data and because of inadequate theoret- 
ical understanding of the coefficients of viscosity for 
these complex mixtures, it was decided to develop 
experimental procedures for measuring their coef- 
ficients of viscosity. It was desired that the tech- 
nique should be simple and accurate for direct mea- 
surement of coefficients of viscosity for all types of 
gases. It should require only small quantities of the 
test gas and should be easily adaptable for relative 
measurements. 

There are several techniques for measuring gas- 
eous viscosity. The prominent ones are summarized 
in the sections which follow. 

General Capillary Method 

In the case of a compressible fluid flowing steadily 
through a perfectly cylindrical capillary tube, the 
volumetric flow rate at the inlet, is given by the 


following expression (refs. 2 to 4): 
v = Krjipi-po) 

8 r}£ c 


where 


r c 

root-mean-square radius of tube 

tc 

length of tube 

n 

coefficient of viscosity of fluid 

Pi 

pressure at inlet 

Po 

pressure at outlet 

This is Poiseuille’s formula for the flow of a gas 
through a tube. (See ref. 5.) Since all the quantities 
in the expression are easily measurable, it provides a 
convenient method for measuring 77 experimentally. 
Equation (1) can be simplified as follows: 


V = A— (2) 

V 

where 

OO 

II 

A more accurate form of this equation would be 
(refs. 6 and 7): 


v = n(i% (1 + ^- B (3) 

where 


8 

capillary bore nonuniformity correc- 
tion factor 

M 

b 

wall-slip correction factor, X 

B 

flow nonlaminarity correction factor 


For very low flow rates (corresponding to Reynolds 
numbers from 100 to 2000) sustained turbulence 
cannot be expected. Under these conditions, B is 
almost zero, and the slip correction is negligible for 
all gases. Hence, equation (2) provides an excellent 
basis for direct determination of viscosity as follows: 

( ^( Air ) A = ( Ap(Air) \ / rj (Test) \ 
Vu(Test) / \Ap(Test)/ \ 77 (Air) / K) 

i/(Test) = „(Air) (j^) (^g*) (5) 

By measuring the differential pressure across the 
same capillary tube for low-volume flow rates 
(< 50 cm 3 /min) for air and the test gas, the viscos- 
ity of the test gas can easily be determined. The 



accuracy of the viscosity value is limited only by the 
accuracies with which F(Air), F(Test), and Ap can 
be measured. 

Oscillating-Disc Method 

This method involves measurement of the period 
and amplitude of the damped oscillations of a suit- 
able solid body. The solid body is suspended from an 
elastic wire first in the test fluid and then in vacuum. 
A simplified expression for the fluid viscosity may be 
written as follows (refs. 8 to 10): 



requirement that no eddy motion develop in the fluid 
flow between the two cylinders. 

The Falling-Body Method 

The basis of the method lies in Stokes’ law 
(ref. 13), which states that the viscous drag W for 
a sphere falling with a uniform velocity in a fluid is 
given by 

W — 67 xr}r s V'p (8) 

where 

r s radius of falling sphere 

Vx terminal uniform velocity 


where 

A, A 0 = damping coefficients in test fluid and 
vacuum, respectively 

T, T 0 = period of oscillation in test fluid and 
vacuum, respectively 


By measuring the terminal velocity of a sphere of 
known dimensions, the viscosity of the test medium 
can be obtained from the following expression: 

” = l gr > (9) 


C = ^ d ^ +d . 2 ) ’ an ex P er i men f a l constant 

Id = moment of inertia of oscillating disc 

r fjj = radius of oscillating disc 

d\,d 2 = distances between oscillating and fixed 
discs 


However, the velocity field around the oscillating 
body is quite complex. This has been the major 
limitation of this method until recently. 


Rotating-Cylinder Method 


In its most widely used variant, the angular de- 
flection 9 of the inner cylinder is noted when the outer 
cylinder is rotated with a constant angular velocity 
oj. If the end effects can be ignored, the coefficient 
of the fluid viscosity is given by (see refs. 11 and 
12 ): 


-KOIcjrl - rf) 
r^rfT 2 uj£i 


(7) 


where 


I c moment of inertia of inner cylinder 

£i length of inner cylinder 

T natural period of vibration 

Tj-, r 0 radii of inner and outer cylinders, 
respectively 

The major limitations of this method are the com- 
plexity of calculation of end-effect corrections and the 


where 


Ps 

density of falling sphere 

Pf 

density of test fluid 

9 

acceleration due to gravity 


The major limitations of this method include the 
problems of frictional charge and the establishment 
of a truly uniform fall through the test fluid. 

From these various techniques, it was decided 
to use a variant of the capillary flow method for 
measuring the viscosities of the test gases. This 
method is particularly suited for complex mixtures, 
which can often be costly and are available only in 
small quantities. 

As indicated previously, the capillary tube vis- 
cometer is based on the application of Poiseuill’s law 
for the steady flow of Newtonian fluids through a cap- 
illary tube. Clearly, all the assumptions made in the 
derivation of the laminar-flow solution cannot be ful- 
filled in its practical design and operation. However, 
a detailed study of all the sources of errors indicates 
that a practical capillary tube viscometer, capable of 
an accuracy of the order of 1 percent, can be devel- 
oped if the following conditions are met: 

1. A laminar steady-state flow condition is 
achieved before the differential pressure and volumet- 
ric flow rates are measured. 

2. The fluid flow rates are kept low (i.e., Reynolds 
number is kept below 2 x 10 3 ). 

3. The effects of compressibility of actual fluids 
flowing through the capillary tube are minimized. 
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(This is possible if the flow rates are kept below 
50 cm 3 /min and if 2a ~ l.) 

The experimental procedures and the results for a 
number of binary and multicomponent mixtures are 
discussed in the following sections. 

Symbols 

B flow nonlaminarity constant 

d distance between oscillating disc and 

fixed disc 

fi mole fraction of zth specie 

I c moment of inertia of inner cylinder 

Id moment of inertia of oscillating disc 

I c length of capillary tube 

li length of inner cylinder 

Mi molecular weight of zth specie 

n number of components in multicompo- 

nent mixture 

p pressure 

Pi pressure at inlet of capillary tube 

p 0 pressure at outlet of capillary tube 

r c root-mean-square radius of capillary 

tube 

rd radius of oscillating disc 

rj radius of inner cylinder 

r a radius of outer cylinder 

r s radius of falling sphere 

T period of oscillation in test fluid 

T a period of oscillation in vacuum 

V volumetric flow rate 

Vx terminal uniform velocity 

W viscous drag 

X wall-slip correction factor 

A p differential pressure (pi — p 0 ) 

Xpc pressure drop across series capillary 
tube 

A pp pressure drop across capillary section 

of flowmeter 

6 capillary bore nonuniformity correc- 

tion factor 


r) coefficient of viscosity 

rji coefficient of viscosity of zth 

component 

6 angular deflection of inner cylinder 

A damping coefficient in test fluid 

A c damping coefficient in vacuum 

p f density of test fluid 

p s density of falling sphere 

4>ij Chapman-Enskog parameter for 

species i and j 

oj angular velocity 

Theory 

The theory of viscosity of gases has been discussed 
extensively by a number of authors. (See refs. 14 
to 17.) In this report, the computational procedure 
used for calculating the viscosity of multicomponent 
gas mixtures is briefly reviewed. 

According to the Chapman-Enskog (C-E) kinetic 
theory of gases (ref. 15), the viscosity of a low- 
pressure multicomponent mixture can be approxi- 
mated by the following expression: 

r](Mixture) = J2 , (10) 

i fjtij 

where 

n number of components 

fi mole fraction of fth specie 

There are a number of methods for estimating the 
parameter <f>ij. We have used Wilke’s approximation 
(refs. 18 to 20) for <pij as follows: 



where 

fi mole fraction of zth component 

Mi molecular weight of ith component 

For nonpolar mixtures, Wilke’s approximation is 
preferable over Brokaw’s treatment (ref. 21) because 


3 



of its simplicity. It has been reported to be accurate 
to within 1 percent for nonpolar binary mixtures. 
For example, for a binary mixture, 1 one obtains the 
following: 


rj (Mixture) = 


Vlfl , mh 
fi + h<t>i2 h + /i02i 


(13) 


It is apparent that the viscosity of the mixture 
is not equal to the weighted mean of the values 
of the viscosity of its components. Equation (10) 
has been used for calculating t) (Mixture) for var- 
ious types of mixtures. The viscosities of pure 
gases have been taken from references 6 and 23. 
The calculational program is described in the next 
section. 


Computational Procedure 

The computer program WILKE is written in 
FORTRAN Version 5 language for the Control 
Data CYBER 170 series digital computer system 
with Network Operating System (NOS) 2.3. The 
program requires approximately 26000 octal lo- 
cations of core storage. A typical run to ana- 
lyze 30 cases, each with 5 component gases, re- 
quires approximately 1 CPU second on the 
CYBER 173. 

Program execution is performed by case, and the 
user may provide data for as many cases as desired. 
For each case, data are input to the program in 
NAMELIST format using the name VAR. Four vari- 
ables must be defined for each case. The integer 
case number (NCASE) is used only for bookkeep- 
ing. The integer number of gases to be considered 
(NGAS) may range from 2 to 12. An integer ar- 
ray (IPOINT) is used for selecting which gases are 
to be included. A value must be provided for each 
of the NGAS gases. Gases are selected by specify- 
ing the appropriate indices from table A. CONC is 
a real array defining the concentration of each gas. 
The sum of all concentration fractions should equal 1, 
and the position within the array for each concentra- 
tion should also follow the gas/index values provided 
above. As an example, a run to analyze two cases is 
presented below. The first case consists of a mixture 
of 30% CH 4 , 50% N 2 , and 20% C0 2 . The second 
case consists of a mixture of 35% CH4, 40% C2H6, 
and 25% C 4 Hi 0 . 

1 For multicomponent mixtures, Buddenberg and Wilke 
(ref. 22) recognized the necessity of introducing the diffusion 
coefficients into the formula for r] (Mixture) to take into 
account the collisions between unlike pairs. Even though this 
modified expression has been reported to predict experimental 
results quite well, it is not always usable because of lack of 
data about diffusion coefficients. 


Gas 

Index 

ch 4 

1 

c 2 h 6 

2 

c 3 h 8 

3 

n 2 

4 

c 4 h 10 

5 

co 2 

6 

0 2 

7 

nh 3 

8 

Ne 

9 

Ar 

10 

H 2 

11 

He 

12 


Table A 


$VAR NCASE = l.NGAS = 3, IPOINT = 1,4,6, CONC = 
0.30, 0.00, 0.00, 0.50, 0.00, 0.20$ 

$VAR NCASE = 2, NGAS = 3, IPOINT = 1,2,5, CONC = 
0.35,0.40,0.00, 0.00,0.25$ 


Program output for each case echoes the case 
number and concentrations of selected component 
gases and gives the viscosity of the mixture. A listing 
of the program appears on the following page. 

Experimental Procedure and Results 

It was the purpose of this study to develop a tech- 
nique for measuring the viscosity of the unknown 
test gas without knowing its physical characteristics. 
This required the use of a dry test meter in series 
with an appropriate mass flowmeter for determining 
the calibration constant of the flowmeter for the test 
gas. Using this calibration constant, the volumet- 
ric flow rate of the test gas through any other mass 
flowmeter based on the same principle can easily be 
determined. As indicated previously, a variant of 
the capillary tube method has been selected for mea- 
suring the viscosity of various types of gaseous mix- 
tures. This method is simple, sensitive, and adapt- 
able to absolute or relative viscosity measurements 
for low-pressure gases. Figure 1 is a schematic of the 
experimental system. The capillary section of the 
flowmeter is made up of two narrow (0.030-cm and 
0.045-cm i.d.) capillary tubes in parallel. The flow 
through each element of the capillary section remains 
laminar under low flow rates. Figure 2 shows the cap- 
illary section of the flowmeter selected for this study. 
The pressure difference accross the capillary section 
of the selected flowmeter is measured for air and then 
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PROGRAM W ILKE! INPUT# OUTPUT# TAPE5 -I NPUT#T A PE6 -OUTPUT) 

DIMENSION MASS(15),C0NC (15)# ETA (15) 

DIMENSION DET A ( 15 ) 

REAL MASS 
REAL K# MU 
DIMENSION S YM B ( 15 ) 

DIMENSION IP0INTU5) 

CHARACTER*10 SYMB 

NAMELIST/VAR/IPOINT#NGAS»CONC#NCASE 

DATA MASS/ 16. #30. #44. #28. #5 8. #44., 32. #17. #20. #40. #2., 4. #3*0./ 
DATA ETA/109.8E— 6#91.0E-6# 80 . OE-6# 175 . 69E-6# 74. 8E-6# 

1 146.63E-6#203.31E-6#97.4E-6#313.81E-6# 

2 222.86E-6#88.73E-6,196.14E-6,3*0./ 

DATA DETA/ 0.330# 0.277# 0.220# 0.454# 0.237# 

1 0. 450# 0. 616# 0. 425# 0.697#0. 704# 0.200# 0.464# 

2 3*0./ 

DATA SYMB/'CH4»#'C2H6«# »C3H8'#'N2«# • C4H10 • # • C02 • # • 02 • # 

1 'NH3'#'NE'»'AR'#»H2'#'HE(4)'#3*' ’/ 

CHANGE THIS CARD FOR A DIFFERENT TEMPERATURE 

DATA T/300./ 

1 FORMAT ( 1H1 ) 

2 FORMAT! ////14X# 'CASE NUMBER '#12///) 

3 FORMAT (///9X» 'COMPONENT '#6X# ‘FRACTION'//) 

4 F0RMAT(9X#A10#6X#F7.4/) 

5 FORMAT! ///3X# 'VISCOSITY OF MIXTURE - »#1PE11.3) 

DO 8 I JK-1# 15 

DETA(IJK) - DETA(IJK)*l.E-6 
8 CONTINUE 
10 CONTINUE 

READ!5#VAR#END-20) 

20 IF!EOF(5).NE.O) GO TO 999 
SUM - 0. 

DO 40 I "1# NGAS 
ANS « 0. 

DO 30 J "1# NGAS 
IFII.EQ.J) GO TO 30 
IND1 - I POINT ( I ) 

I ND2 - IP0INT1J) 

ETA1 - ETAUND1) + DETA!IND1)*(T - 293.15) 

ETA2 - ETA ( I ND2 ) ♦ D ET A { IND2 ) * ( T - 293.15) 

PHI - 1. ♦ SORT !ETA1/ETA2)*(MASS ( IND2) /MASS! IND1) ) 

1 **0.25 
PHI ■ RHI*PHI 

PHI - PHI/SQRT! 8. *11. + MASS ! IND1 ) /MASS ( IND2 )) ) 

ANS - ANS + CONC ( IND2 ) *PHI 
30 CONTINUE 

SUM « SUM + ETA1/(1. + ANS /C ONC ! IND1 ) ) 

40 CONTINUE 
WRITE ( 6# 1 ) 

WP I TE ( 6# 2 ) NCASE 
WR I T E ( 6# 3 ) 

DO 50 I -1# NGAS 
IND - I POINT! I ) 

WPITE(6#4) SYMB(IND)#CONC( IND) 

50 CONTINUE 

WRITE (6# 5 ) SUM 
GO TO 10 
999 CONTINUE 
STOP 
END 
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for the test gas. 2 Care is taken to keep the flow rates 
for the air and the test gas below 50 cm 3 /min to en- 
sure laminarity of the gas flows through the capillary 
section of the flowmeter. Figure 3 shows a correlation 
between the pressure difference across the test capil- 
lary section and the flow rate through it. It is obvi- 
ous that, as would be expected for laminar flows from 
equation (2), the relationship between the flow rates 
and the corresponding differential pressures is linear 
for flow rates from 10 to 50 cm 3 /min. In this experi- 
mental arrangement, the flowmeter provides data for 
determining the volume flow rates, as well as a well- 
characterized capillary tube section for differential 
pressure measurement across it. It is this dual role 
of the flowmeter that constitutes the uniqueness of 
the method developed in this study. The results for 
binary and more complex mixtures are summarized 
below. All measurements were made at room tem- 
perature and at a pressure of 1 atmosphere. Figure 4 
is a photograph of the viscosity measurement station. 

Purity of gases. Table B lists the components and 
their purity as quoted by the suppliers. 


Gas 

Purity 

Air 

Ar 

C0 2 

n 2 

ch 4 

c 2 h 6 

c 3 h 8 

c 4 h 10 

Dry grade; Dew point = — 80°F 
99.998% 

99.50% 

99.998% 

99.9995% 

99.00% 

99.95% 

99.50% 


Table B 


Preparation of binary mixtures. The binary mix- 
tures were prepared by mixing appropriate volume 
flow rates of the individual gases. This procedure 
ensured complete interdiffusion of the two compo- 
nents. Selected gases were mixed in various ratios to 
produce a total flow rate from 20 to 50 cm 3 /min. Dif- 
ferential pressure measurements across the test cap- 
illary element of the gas flowmeter were made with 
a differential capacitance gage. Viscosity values for 
the test mixtures were then calculated by using equa- 
tion (5). The results are summarized in table I. 


2 Pressure differences were measured across flowmeter 1 

and across series capillary tube section C-l (fig. 1). The 
values of the viscosity of pure gases calculated on the basis 
of these two independent pressure differences were found to 
be equal within approximately ±1 percent. It was therefore 
decided to eliminate the series capillary section from all 
further measurements. 


Preparation of complex mixtures. The more com- 
plex mixtures were prepared by a commercial vendor 
according to our specifications. The vendor prepared 
the mixtures by gravimetrically mixing different mass 
fractions, which were later expressed in appropriate 
mole fractions. The test mixtures were kept in their 
containers for several weeks before use to ensure com- 
plete mixing. 

Four different mixtures containing different vol- 
ume fractions of C x H y hydrocarbons as well as N 2 
and CO 2 were prepared. Their viscosities were mea- 
sured using the same procedure described for binary 
mixtures in the preceding section. The results are 
summarized in table II. 

Natural gas. A test sample of the natural gas 
was obtained from the Virginia Natural Gas Com- 
pany. Its coefficient of viscosity was measured by 
using the same procedure as for the binary and com- 
plex mixtures. The test sample was also analyzed for 
its composition. The composition data are summa- 
rized in table III. These values are averages of two 
measurements of the same sample. The accuracies of 
individual gas chromatographic analyses are of the 
order of 1 percent (relative). Using these data, the 
coefficient of viscosity of the sample was calculated, 
and the results are summarized in table IV. 

Pure gases. Even though the main purpose of 
this study was to develop a technique for measuring 
the viscosity of multicomponent gas mixtures, it was 
deemed desirable to measure the viscosity of selected 
pure gases as a test of the accuracy of the proposed 
technique. These selected pure gases were treated as 
unknown test gases, and their viscosities were mea- 
sured by using the same procedure as was used for 
the multicomponent mixtures. The present results, 
along with the values reported in the literature are 
summarized in table V. It is apparent that the two 
sets of data are in excellent agreement. All measured 
values are strongly dependent on the viscosity of air 
at room temperature and pressure. The viscosity of 
dry air has been taken to be 185.6 p P. (See refs. 6 
and 23.) 

Discussion 

As seen from the data in table I, the agreement 
between the measured and calculated values of vis- 
cosity for binary mixtures is quite good. Also in- 
cluded in the table are some of the values reported by 
other authors. These values are also in good agree- 
ment with experimental values. Table II summa- 
rizes the experimental and calculated values for four 
complex mixtures prepared to simulate U.S. natural- 
gas samples. In these cases, the agreement is less 
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satisfactory — possibly because of uncertainties in the 
mixture compositions. To further resolve this is- 
sue, the viscosity of a Virginia Natural Gas sam- 
ple was measured. This sample was also analyzed 
for its composition. The measured value of its vis- 
cosity using the present technique (independent of 
gas composition) is 106.7 fiP (±2.0): the expected 
value was 107.9 fiP (±2.0). These values should be 
compared with a calculated value of 111.6 /iP (±2.0) 
based on the listed gas composition. These results 
indicate that ^(calculated) values for complex mix- 
tures containing hydrocarbons are slightly different 
from ^(Measured) values, possibly because of the 
limitations of Wilke’s approximation of the classi- 
cal dilute-gas state model. A direct measurement 
of the viscosity of the unknown test mixture may 
therefore be preferable to the one computed on the 
basis of its measured composition, since the limita- 
tions of Wilke’s approximation of the dilute-gas state 
model make the latter less accurate for complex gas 
mixtures. 

Concluding Remarks 

A variant of the general capillary method for mea- 
suring the viscosity of unknown gases has been de- 
veloped. It is based on the use of a thermal-mass 
flowmeter whose critical sensing element simultane- 
ously provides a capillary section for direct measure- 
ment of pressure drop. It is this dual role of the 
thermal-mass flowmeter that constitutes the unique- 
ness of the new technique. It does not require any 
prior information about the composition or other 
physical properties of the test gases if a dry test meter 
is used for initial calibration constant determination. 
It is particularly suited to very complex hydrocarbon 
mixtures where classical theory limitations as well 
as compositional errors make theoretical calculations 
less reliable. The results for pure gases indicate that 
the technique is accurate to within ±1 percent. In 
fact, the only limiting factors appear to be the accu- 
racies with which volume flow rates and differential 
pressure can be measured. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
February 20, 1986 
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TABLE I. SUMMARY OF COEFFICIENTS OF VISCOSITY FOR BINARY MIXTURES 

[T = 300 K\p = 1 atm] 


Gas mixture, 
percent volume 

Volume 
flow rate, 
cm 3 /min 
(±1%) 

Differential 

pressure, 

torr 

(±0.25%) 

Coef 

icient of viscosity, pP 

^(Literature) 

T) (Measured) 3 


Air 

48.9 

0.693 

185.6 



24.9% N 2 + 75.1% C0 2 

46.6 

.555 

c 158.2 ± 1.6 

156.0 

155.6 

50.2% N 2 + 49.8% C0 2 

45.8 

.573 

c 165.7± 1.7 

163.8 

162.3 

75.0% N 2 + 25.0% C0 2 

46.1 

.597 

c 172.9± 1.7 

169.6 

169.8 

Air 

48.9 

0.708 

185.6 



25.1% N 2 + 74.9% Ar 

39.6 

.683 

d 215.8±0.2 

221.1 

216.9 

49.8% N 2 + 50.2% Ar 

39.0 

.623 

d 204.3 ± 0.2 

204.9 

204.9 

75.3% N 2 + 24.7% Ar 

38.3 

.560 

d 191.5±0.2 

187.3 

192.1 

87.7% N 2 + 12.3% Ar 

35.9 

.520 

d 186.1±0.2 

185.5 

185.2 

Air 

45.9 

0.689 

185.6 



24.4% Ar + 75.6% C0 2 

51.4 

.690 

e 167.5± 1.7 

166.1 

164.8 

48.8% Ar + 51.2% C0 2 

51.0 

.762 

e 186.0± 1.9 

184.9 

181.9 

74.1% Ar + 25.9% C0 2 

50.6 

.848 

e 205.5 ± 2.1 

207.5 

202.8 

Air 

48.8 

0.671 

185.6 



25.9% N 2 + 74.1% CH 4 

44.4 

.435 


132.2 

130.8 

50.8% N 2 + 49.2% CH 4 

44.1 

.486 


148.8 

147.8 

75.1% N 2 + 24.9% CH 4 

45.6 

.554 


164.0 

162.9 

Air 

49.6 

0.708 

185.6 



25.9% N 2 + 74.1% C 2 H 6 

45.2 

.392 


112.1 

109.4 

50.1% N 2 + 49.9% C 2 H 6 

45.7 

.448 


126.7 

127.4 

75.4% N 2 + 24.6% C 2 H 6 

49.3 

.566 


148.5 

149.9 

Air 

49.6 

0.708 

185.6 



25.0% C0 2 + 75.0% CH 4 

44.2 

.435 

e 126.5± 1.3 

128.0 

127.6 

50.0% C0 2 + 50.0% CH 4 

44.6 

.473 

e 138. 1 ± 1.4 

137.9 

138.0 

75.0% C0 2 + 25.0% CH 4 

46.1 

.515 

e 145.4± 1.5 

145.3 

145.1 

Air 

49.6 

0.708 

185.6 

109.3 


25.0% C0 2 + 75.0% C 2 H 6 

43.9 

.369 



108.1 

50.0% C0 2 + 50.0% C 2 H 6 

45.4 

.436 


124.9 

121.9 

75.0% C0 2 + 25.0% C 2 H 6 

48.4 

.503 


135.1 

135.9 

Air 

49.6 

0.708 

185.6 



24.8% C0 2 + 75.2% C 3 H 8 

36.2 

.276 


98.7 

95.0 

49.9% C0 2 + 50.1% C 3 H 8 

43.9 

.387 


113.9 

109.9 

75.0% C0 2 + 25.0% C 3 H 8 

45.5 

.462 


131.3 

127.9 


“See footnotes on page 10. 
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TABLE I. Concluded 
[T = 300 K-,p = 1 atm] 


Gas mixture, 
percent volume 

Volume 
flow rate, 
cm 3 /min 
(±1%) 

Differential 

pressure, 

torr 

(±0.25%) 

Coefficient of viscosity, //P 

rj (Literature) 

//(Measured) 0 

//(Calculated) 6 

Air 

48.9 

0.702 

185.6 



25.0% CH 4 ± 75.0% C 2 H 6 

39.1 

.298 


98.5 

98.0 

50.0% CH 4 + 50.0% C 2 H 6 

38.9 

.307 


102.1 

102.0 

74.9% CH 4 + 25.1% C 2 H 6 

40.2 

.330 


106.1 

106.6 

Air 

48.9 

0.702 

185.6 



25.8% CH 4 + 74.2% C 3 H 8 

33.8 

.224 

f 86.8 9 86.0 

85.7 

85.1 

50.0% CH 4 + 50.0% C 3 H 8 

39.8 

.285 

*92.9 9 93.0 

92.7 

93.5 

74.9% CH 4 ± 25.1% C 3 H 8 

39.2 

.302 

o 

CN 

o 

/-H 

b- 

o 

o 

T— 1 

99.5 

101.4 

Air 

48.9 

0.702 

185.6 



25.0% C 2 H 6 ± 75.0% C 3 H 8 

24.7 

.165 


86.4 

85.0 

51.0% C 2 H 6 + 49.0% C 3 H 8 

36.8 

.254 


89.2 

88.0 

74.9% C 2 H 6 + 25.1% C 3 H 8 

39.8 

.283 


91.9 

91.0 


“The uncertainties in //(Measured) values, determined by the accuracies of the gas flow rate and 
differential pressure measurements, are within the range of ±2.5 percent in the present study. 

6 The uncertainties in //(Calculated) values, stemming from uncertainties in the mixture composi- 
tion, are within the range of ±2 percent. 

“These values are the results of measurements of absolute viscosity using an oscillating disc 
viscometer (ref. 22). 

d These values are the results of relative measurements of viscosity using an oscillating disc 
viscometer (ref. 23). 

“These binary mixture viscosity measurements were made using an oscillating disc viscometer 
(ref. 24). 

•^These data are cited in reference 15 for comparison between the experimental data and the 
calculated values. They are in agreement with Wilke’s approximation to an accuracy of ±1 percent. 

9 These are absolute values obtained using a high-pressure capillary tube viscometer. The repro- 
ducibility of the data has been reported at ±0.25 percent (ref. 25). 
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TABLE II. SUMMARY OF COEFFICIENTS OF VISCOSITY FOR COMPLEX MIXTURES 


r?( Air) = 185.6 pP; 
T = 300 K\ 
p = 1 atm 


Gas mixture, 
percent volume 

Volume 
flow rate, 

Differential 

pressure, 

torr 

Coefficient of viscosity, pP 

(± 1 % relative) 

cm 3 /min 

(±0.25%) 

7 ? (Measured ) 0 

? 7 (Calculated) b 

Air 

44.6% C 2 H 6 + 25.9% C 3 H 8 
+ 10 . 8 % C 4 H 10 + 

5.1% C0 2 + 13.6% N 2 

48.9 ±0.5 
31.8 ± 1.0 

0.682 

.220 

92.0 ± 1.5 

97.0 ±2.7 

53.0% CH 4 + 26.1% C 2 H 6 
+ 18.4% C 3 H 8 + 

1.4% C0 2 + 1.0% N 2 

39.3 ± 1.2 

.318 

105.3 ± 1.1 

100.7 ±2.8 

50.8% CH 4 + 17.1% C 2 H 6 
± 24.9% C 3 H 8 + 

3.3% C0 2 + 3.8% N 2 

36.8 ± 1.4 

.298 

105.7 ±0.7 

102.1 ±2.9 

51.5% CH 4 ± 21.2% C 2 H 6 
+ 23.8% C 3 H 8 + 

1.7% C0 2 + 1.7% N 2 

44.6 ± 0.9 

.308 

91.7 ±0.7 

100.0 ± 2.8 


“The errors cited on ^(Measured) are standard deviations on 5 measurements for each mixture. 
They do not include the effects of uncertainties in flow rates and differential pressure measurements. 

^The uncertainties cited on ^(Calculated) represent the effects of uncertainties in the mixture 
composition. 
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TABLE III. COMPOSITION OF NATURAL 
GAS SAMPLE 

[Ref. 21] 


Component gas 

Volume percentage 

ch 4 

96.52 ±0.20 

c 2 h 6 

2.32 ± 0.08 

c 3 h 8 

0.40 ± 0.03 

C 4 H 10 

0.14 ±0.02 

C 5 H 12 

< 0.05 

C 6 H 14 

< 0.05 

C0 2 

0.22 ±0.02 

n 2 

0.39 ±0.06 


TABLE IV. SUMMARY OF COEFFICIENTS OF VISCOSITY FOR 

NATURAL GAS 


77 (Measured), /xP 

^(Calculated), 0 ^P 

Based on gas composition 
summarized in table III b 

Independent of 
gas composition 0 

Based on gas composition 
summarized in 
table III 

107.9 ± 2.0 

106.7 ± 2.0 

111.6 ±2.0 


“The uncertainty cited on ^(Calculated) represents the effects of uncertain- 
ties in the mixture composition. 

b The value measured on the basis of known gas composition is based on the 
flow-rate calibration constant calculated for that composition. The error cited 
on it is the standard deviation based on 6 measurements. 

“The measured value listed in the second column has been determined using 
the experimentally determined flowmeter calibration constant for the test gas. 
The error cited on it is the standard deviation based on 6 measurements. 





TABLE V. SUMMARY OF COEFFICIENTS OF VISCOSITY FOR PURE GASES 


[T = 300 K\ p = 1 atm] 


Gas 

Volume 
flow rate 
cm 3 /min, 
(± 1 %) 

Differential 

pressure 

torr 

(±0.25%) 

Coefficient of viscosity, fiP 

Measured , 0 

(average) 

Reported 6 (refs. 6 and 23) 
(± 2 %) 

Air 

48.9 

0.708 


185.6 

n 2 

32.3 

.449 

178.1 ±0.5 

178.6 

C0 2 

35.6 

.412 

150.3 ±0.1 

149.9 

Ar 

41.0 

.723 

228.2 ± 0.2 

227.6 

ch 4 

30.9 

.260 

112.7 ±2.2 

112.0 

c 2 h 6 

37.1 

.266 

95.7 ±2.2 

94.5 

c 3 h 8 

30.1 

.198 

83.9 ± 1.6 

82.5 

C 4 H 10 (iso) 

23.5 

.141 

76.6 ± 1.0 

75.7 


“The errors cited on the measured values are standard deviations based on 
3-5 measurements. 

b The uncertainties on the reported values fall within the range of ±2 percent. 




Air 

Test gas 1 
Test gas 2 



Atmosphere 

Atmosphere 

Atmosphere 


Figure 1. Schematic of experimental system for measuring viscosities of atmospheric-pressure gases. 


Meter 



Figure 2. Thermal flowmeter with laminar flow element. 
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Figure 3. Flow rate versus differential pressure across test capillary section of flowmeter. 
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Figure 4. Photograph of experimental viscosity measurement station. 
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